Motion of a falling object 
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Abstract 



fi , A simple setup was assembled to study the motion of an object while it falls. The setup was 

used to determine the instantaneous velocity, terminal velocity and acceleration due to gravity. 
^>^, Also, since the whole project was done within $20 it can easily be popularized. 

1 Introduction 

V 

\^ • One of the first topic students of physics learn is the concept of bodies falhng towards the center 

^^ I of the earth with an acceleration (g) . This acceleration due to gravity is usually determined in 

undergraduate laboratories by measuring the oscillation time period of a simple pendulum. Some 
students do wonder as to why acceleration due to gravity is not just determined by measuring 

^^ ■ the time taken by a freely falling body to fall through a given height. Trivial calculation using 

f^ I the equation 

^ '. 1 

O, s^ut+-gt^ (1) 

c/2 ; 

j^. shows that the body has to drop across an appreciable height to increase the body's time of flight. 

'" ^ ' This is required to minimize the error in measurement of the fall's time. Historically this would 

explain why Galileo did his experiment from the Tower of Pisa. Over the last three decades 
various contraptions have been designed to measure the time of flight in laboratory conditions 

^^ . i.e. drop through limited height H-IH. This has been made possible by the advent of cheap and 

^ ' easily available electronic circuitry capable of measuring accurately in milliseconds. In most of 

these experiments, the body is released by either a lever or electromagnet. The position of the 
release mechanism would be taken as s=0 and the initial velocity on release would obviously be 
zero. A second sensor kept at a distance one meter below the release point (s=0) would stop 
the timing circuit when the body crosses it. The acceleration due to gravity is calculated using 
eqn(^. However, this scheme introduces an error in time 't' due to the uncertainty in the time of 
release. In Lindemuth's Jt^ work, error was introduced due to the body having an initial velocity 



by the time it arrived to the first sensor that initiated the clock. In all schemes discussed, there 
was an error either in the measurement of time of fall or knowledge of initial velocity. Also, 
these experiments relay on the a prior knowledge/ validity of eqn(Q. To overcome this problem 
a method devised in use is to allow a metallic ball to drop in between two parallel wires in which 
high tension current flows in pulses. Sparks generated by the traveling ball is recorded on a 
photographic strip. In another method two test masses to required to determine acceleration 
due to gravity ^. These techniques either require high (unsafe) voltages or the contraptions 
require a certain amount of skill in designing which would not be easily available in standard 
under-graduate colleges of developing countries. 

We have designed a simple setup to measure the acceleration due to gravity. It was designed 
keeping in mind minimum skill should be required to set it up and use it as also the validity 
of eqnQ would not be taken for granted. In the following sections we elaborate our setup 
designed by using LEDs, photodiodes, op amps and microprocessor and comment on the results 
we obtained. 



2 Apparatus 
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A hollow tube (which we shall hence refer to 
as the " drop tower'^ ) of length Imeter housed 
eight LEDs (Light Emitting Diodes) and eight 
LDRs( Light Dependent Resistances). They 
were arranged along the length of the drop tower 
facing each other (see fig 01 . The distance be- 
tween two neighboring LEDs (and in turn LDRs) 
was 10cm. The first pair of LED-LDR was kept 
at a distance 10cm from the top edge of the drop 
tower. This was to make sure that stray light 
from the room did not lead to false triggering 
of the circuitry. 

As the name suggests, LDRs are light light 
sensitive devices and hence are popularly used 
as optical sensors. Their resistance depends on 
the amount of light falling on it. When light 
is incident on the LDR, it's resistance is low. 
Thus, when a falling object cuts the lights path, 
no light falls on the LDR and it's resistance in- 
creases. This variation in resistance when con- 
verted into voltage can be (wave-)shaped into 

pulses which then can be fed to a computer or microprocessor for identifying and storing the 
information as to the instant the object passed the LDR. The required circuitry used for wave- 
shaping is shown in fig©- The LDR resistance is converted into a corresponding voltage using 
a voltage divider circuit. The voltage across the resistance, Ri, is high (Vh) when light falls 
on the LDR since it's resistance diminishes. Conversely, if the LDR is in shadow region it's 
resistance is high and the voltage across Ri goes low (Vl). This voltage is given to the op amp 



Figure 1: Setup used to measure the time of 
fall of a falling object. 



comparator circuit. The variable voltage is set at a value between Vl and Vh- This ensures 
that the output of the LM741 op amp swings to +12v when light falls on the LDR and to -12v 
when shadow falls on it. This voltage is made TTL/microprocessor compatible (vary between 
and +5v) using a 4.7v zcncr diode. 
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Figure 2: Wave-shaping circuit used to convert variation ofLDR's resistance into square pulses. 
The DC supplies used were Vi = 5v and V2 = 12v and the resistances in the circuit were 
Ri = 3.3K51, R2 = lOKfJ and R3 = IKil. Diodes used were the usual 1N4001 rectihcation 
diodes. 

Now we explain the final block of our setup which is a IC8085 microprocessor training kit. 
To follow in detail this section requires an understanding of IC8085 and can be acquired from 
Gaonkar's book. A standard 8085 microprocessor training kit comes with a programmable 
input/output device (IC8155) which has as many as twenty two input pins (divided into Port 
A with 8 pins, Port B with 8 pins and Port C with 6 pins). The output of each wave-shaping 
circuit can be fed to one of these pins. We have only used eight pins of Port A. This was done 
keeping the following in mind 

(i) Optimize the amount of hardware (wave-shaping circuits) , 

(ii) eight data points would be reasonably enough to fit a curve reliability and 



(iii) the circuit then could easily be extended for interfacing with a computer since the com- 
puters parallel port can only input eight bits. 
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Figure 3: Setup used to measure the time of fall of a falling object. 

When no object cuts the hght path, all the pins of Port A (of IC8155) are high and informs 
the microprocessor by sending the word FFh (the subscript 'H' means that it is the hexa repre- 
sentation of binary number 1111 1111). This word changes when a object cuts the light's path. 
For example on cutting the first LED-light's path the word sent to the microprocessor is 1111 
1110 (FEh) and when the object cuts the second LED-light's path the microprocessor receives 
the word 1111 1101 (FDh) and so forth as the object continuous to fall. The microprocessor 
is programmed to run a counter (clock) from the negative edge of one input signal to that of 
the second (see figEl- On receipt of the second input signal's negative edge, the microprocessor 
stores the counter's count in it's memory and restarts the counter. This contributes in an error 
(time measured would be lower than actual) in the measurement of the time of flight between 
the next two pair of LED-LDR. The counts stored in the memory is converted to time in seconds 
using the formula 



ti 
^2 to is 



(30n- 
(30n- 



34) X 0.326/is 
41) X 0.326/is 



Notice both the wave-shaping circuit and the software are designed to minimize any delay in 
response by the cadmium sulphide LDR used. The assembly language instruction set used for 
the project is listed at the end. 

3 Results and Discussions 



The drop tower was kept vertical and attempt was made to make sure that the objects would 
fall in a straight line. A plumb line was used to ascertain the alignment. For conducting this 
experiment we used four different balls whose physical features, dimensions and mass are listed 
in Table 1. 

The balls were dropped from the top edge of the drop tower. We made sure that the ball in 
use would cut all eight light beams. As soon as the ball crosses the first LED, the microprocessor 
starts it's clock. This point would act as the co-ordinates (0,0). The microprocessor stores the 
counter's count as the ball passes the second LED (the count in turn would be used to calculate 
time). This gives the co-ordinates (t2i, 0.1) where the subscript '21' indicates this as the time 



taken for the fall between the first and second LED. Since, the LEDs are 10cm apart, in meters we 
have s=0.1cni. In this manner, eight data points are generated and stored in the microprocessor. 
The experiment was repeated as many as 200 times with each ball. 



Smooth Rubber Ball 




Figure 4: The time versus position graph for BaU 1 (smooth rubber hall) and Ball 2 (ping-pong 
hall). 

FigQ shows two curves each for the results obtained with Ball 1 and Ball 2. As can be seen 
the data lie on a parabola. However, before fitting a second polynomial that would represent 
eq([Pl, one has to remember that eql^l represents a object falling freely under the action of 
gravity. No other forces should be acting on it. However, there are two possible forces acting in 
the opposite direction to gravity when the object is dropped, namely forces of buoyancy and air 
resistance. The contribution due to air resistance is proportional to the square of the velocity 
with which the object is moving. The general equation of motion in this case is given as 

<Ps 
'm—jj = (m - m.air)g ~ kv"^ 

where niair is the mass of the air displaced by the ball and 'k' is the proportionality constant. 
Since, the dimensions of the balls are small, the volume and in turn the mass of air displaced is 
very small. Hence, contributions of the buoyancy force can be neglected. Hence, the equation 
of motion reduces to 



di2 



mg 



kv^ 



(2) 



Table 1. 



Ball Name 


Description 


Mass (grams) 


Radius (cm) 


Balll 
Ball 2 
Balis 


Smooth Rubber Ball 

Ping Pong Ball 

Nylon Ball 


19.00 
1.75 
5.50 


1.4 
1.9 
1.9 



Before proceeding, for completeness it should be stated here that the force acting as a resistance 
on the falling body can either be proportional to it's velocity (given by Navier-Strokes equation, 
F = 67r/Ltrv) or the square of the velocity. This depends if motion of the viscous medium around 
the falling body remains laminar or not and is denoted by the Reynolds number (R) 



R 



dpv 



where this number characterises the media (viscosity, /i and density, p) through which the body 
(of diameter, d) falls with velocity, v. In a very viscous media (/i large) , the body would face large 
resistance and would fall slowly (v small). Thus the Reynolds number would be less than one. 
Under such conditions the N A vier- Strokes equation should be considered. However, as in our 
case, the body is falling through (dry) air, whose viscosity is small (around 1.9 x 10~^kg/m.sec) 
with density 1.225kg/m'^. The Reynolds number works out be « 3000, i.e. R ::^ 1. The drag 
force acting on the body is then given as 

F = ^pACdv^ 
= kv^ 



or k = ^C^pA, where 'A' is the area of cross-section and Cd is the drag coefficient. The drag 
coefficient depends on the shape, surface characteristics and the Reynolds number JO]. The 
value of Cd for smooth spheres as a function of Reynolds number has been listed in Landau 
and Lifshitz book ^Jj. For R=3000, Cd ~ 0.5. For the radii of the balls used in our experiment 
(see Table 1), k works out to be « 3.8 x 10~^(Ns^/m^). However, we shall not use this value 
immediately but at the end compare our experimentally determined values with this theoretical 
value. 

The solution of eqn(|2Jl and it's power approximation has been detailed by Lindemuth and is 
given as 
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(3) 



where 



Vt 



is the terminal velocity the falling body would achieve when it moves with no acceleration, i.e. 
net force acting on it is zero in the case of free fall else when all other higher terms cancel out. 
The coefficient of t^^ represents the rate of change in acceleration (dg/dt). 

Thus, under the conditions in which this experiment was done, the data points were fit with 
a third order polynomial (eqn|2l. Curves were fit to these data points using standard and freely 
available software "Curxpt v3.0". 



Table 2. 



Ball Name 




ui (m/s) 


g (m/s^) 


dg/dt (m/s^) 


co-relation 


Balll 


I 


1.559 


9.7998 


8.0e-9 


0.998 


II 


1.559 


9.7998 


4e-9 


0.998 


Ball 2 


I 


1.363 


9.798 


0.643 


0.999 


II 


1.368 


9.7745 


0.731 


0.999 


Ball 3 


I 


1.549 


9.7998 


0.418 


0.998 


II 


1.445 


9.7822 


0.268 


0.999 



All the curves for the 200 odd experiments done with each ball lie within the two curves 
shown in figQ). The coefficient of t, t^ and t'^ as returned by the software was used to determine 
the ball's velocity (ui) as it crosses the first sensor where the ball has traveled through lO+cm 
(10cm from the drop tower's top edge to the first sensor and the unknown height from the edge 
from where the ball was dropped) and it's acceleration (g). Results for the graphs shown in 
figlQl) are listed in Table 2. The error in 'g' is 1 in part of 10'^ with as small as +7/isec (max) 
error in the measurement of time of flight. The acceleration due to gravity in Delhi can be 
calculated using the equation 



9 



9.780327[1 + 0.0053024sm2L - 0.000058sm^2L] 



(4) 



where 0.499455 (28°37'N) is the latitude (L) of Delhi in radian. This gives the acceleration due 
to gravity in Delhi as 9.7918m/s^. Our experimental value of 'g' is in good agreement with 
the theoretical value. Notice that the velocity (ui) is different in each case (see Table 2). This 
is expected since even with the best of care taken, the balls would have been dropped from 
different heights (position from the drop tower's edge). 

The power series of eqn© suggests that a knowledge of the coefficients should also enable 
us to estimate the objects terminal velocity. However, the power approximation leads to the 
situation that if the body's initial velocity is zero then the terms associated with terminal velocity 
also collapses. Also, curve fitting the original solution relating distance to time proves tedious 
due to the existence of logaritm and cosine hyperbolic terms. However, the functional relation 
between velocity and time for a body falling through a viscous medium is given as |4, 
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— t + tanh — 
vt \vt 



V = VTtanh 
When the body's initial velcoity is zero, the above expression reduces to 



V — VTtanh 



-t 



Vt 



(5) 



One can either use the coefficients 'g' and 'dg/dt' in 

V ^ gt~3—t 
at 



Table 3. 



Ball Name 




VT (m/s) 


t (sec) 


s (m) 


k = mg7v:^ xlO-4(NsVm") 


Balll 


I 


35.11 


3.58 


62.9 


1.27 


II 


35.11 


3.58 


62.9 


1.27 


Ball 2 


I 


6.32 


0.58 


1.77 


4.3 


II 


5.89 


0.53 


1.52 


4.92 


Balis 


I 


7.97 


0.74 


2.87 


8.47 


II 


9.97 


0.94 


4.6 


5.4 



to generate data for curve fitting eqn© or find the gradient along the curves shown in figQ). 
Both methods would introduce an error abid small since the fitting as represented by co-relation 
factor (see Table 2) was very good. The terminal velcoities thus evaluated is listed in Table 
3. Before concluding, it would be interesting to compute the time a body would take to attain 
terminal velcoity and the distance through which it would fall to attain this velocity. For this 
one would have to solve the quadratic equation 

, 0^9.2 

and substitute the time taken thus obtained in the equation 

2^ dt 

The data of Table 3 indicates that the terminal velocity, the time in which the object attains 
this velocity and the distance through which it has to travel for attaining the same all strongly 
depends on the objects mass. The drag coefficient also depends on the surface area (and hence 
diameter) hence comparison was sort to be done on balls with similar diameter and different 
mass. Obtaining balls of different mass and similar dimensions was not easy. However, the 
designed experimental setup proved handy in estimating the acceleration due to gravity. Above 
which, it also proved to be a rich source of information on the motion of an object while falling. 
Providing ample scope of and experience in data analysis, the setup gives an idea of the terminal 
velocity, time taken to achieve it and the distance through which the body has to travel before 
attaining it. All this is possible eventhough the length of the drop tower is less than the distance 
through which object has to fall to attain vt- On a final note it must be mentioned that this 
project was done in a graduate college laboratory, with a small budget of $20. 
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Program required. 



Mem. add. 


Hex Code 


Instruction 


ClOO 


31 


LXISP 


ClOl 


00 


00 


C102 


C6 


C6 


C103 


21 


LXIH 


C104 


00 


00 


C105 


CO 


CO 


C106 


3E 


MVI A 


C107 


00 


00 


C108 


D3 


OUT 08 (CR) 


C109 


08 


08 


ClOA 


DB 


IN 09 (PA) 


ClOB 


09 


09 


ClOC 


A6 


ANAM 


ClOD 


C2 


C2 


ClOE 


OA 


OA 


ClOF 


CI 


CI 


Clio 


23 


INXH 


cm 


11 


LXI D (initiate counter) 


C112 


00 


00 


C113 


00 


00 


C114 


13 


INXD 


C115 


DB 


IN 09 (PA) 


C116 


09 


09 


C117 


A6 


ANAM 


C118 


C2 


JNZ 


C119 


11 


11 


CllA 


CI 


CI 


CllB 


23 


PUSH D 


cue 


D5 


INX H 


CUD 


C3 


JMP 


CUE 


11 


11 


CllF 


CI 


CI 
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Data used by main program. 



Mem. add. 


Data 


COOO 


01 


COOl 


02 


C002 


04 


COOS 


08 


C004 


10 


C005 


20 


C006 


30 


COOT 


80 



